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We study both linear and nonlinear surface waves localized at the interface separating a left- 
handed medium (i.e. the medium with both negative dielectric permittivity and negative magnetic 
permeability) and a conventional (or right-handed) dielectric medium. We demonstrate that the 
interface can support both TE- and TM-polarized surface wmes-surface polaritons, and we study 
their properties. We describe the intensity-dependent properties of nonlinear surface waves in three 
different cases, i.e. when both the LH and RH media are nonlinear and when either of the media is 
nonlinear. In the case when both media are nonlinear, we find two types of nonlinear surface waves, 
one with the maximum amplitude at the interface, and the other one with two humps. In the case 
when one medium is nonlinear, only one type of surface wave exists, which has the maximum electric 
field at the interface, unlike waves in right-handed materials where the surface-wave maximum is 
usually shifted into a self-focussing nonlinear medium. We discus the possibility of tuning the wave 
group velocity in both the linear and nonlinear cases, and show that group- velocity dispersion, which 
leads to pulse broadening, can be balanced by the nonlinearity of the media, so resulting in soliton 
propagation. 

PACS numbers: 78.68.+m, 42.65.Tg 



I. INTRODUCTION 

Novel physical effects in dielectric media with both 
negative permittivity and negative permeability were 
first analyzed theoretically by Veselago [1] who predicted 
a number of unusual phenomena including, for exam- 
ple, negative refraction of waves. Such media are usually 
known as left-handed (LH) media since the electric and 
magnetic fields form a left handed set of vectors with the 
wave vector. The physical realization of such LH me- 
dia was demonstrated only recently [2] for a novel class 
of engineered composite materials, now called LH meta- 
materials. Such LH materials have attracted attention 
not only due to their recent experimental realization and 
a number of unusual properties observed in experiment, 
but also due to the expanding debates on the use of a 
slab of a LH metamaterial as a perfect lens for focusing 
both propagating and evanescent waves [3]. 

The concept of a perfect lens was first introduced by 
Pendry [4] , who suggested the idea that a slab of a loss- 
less negative-refraction material can be used for creating 
a perfect image of a point source. Although the concept 
of a perfect lens is a result of an ideal theoretical model 
employed in the analysis [4], the resolution limit of a 
LH slab was shown to be independent of the wavelength 
of the electromagnetic wave (but can be determined by 
other factors including losses, spatial dispersion, and oth- 
ers), so that the resolution can be indeed much better 
than the resolution of a conventional lens [5] . 

The improved resolution of a LH slab and the cor- 



responding amplification of evanescent modes even in 
a lossy LH material, can be understood from simple 
physics. Indeed, the near field of an image, which can not 
be focused by a normal lens, can be transferred through 
the slab of a LH material due to the excitation of sur- 
face waves (or surface polaritons) at both interfaces of 
the slab. Therefore, as the first major step in the under- 
standing of the amplified transmission of the evanescent 
waves, as well as other unusual properties of the LH ma- 
terials, it is important to study the properties of different 
types of surface wave that can be excited at the interfaces 
between LH and conventional (or right-handed, RH) me- 
dia. Some preliminary studies in this direction included 
calculation of the linear dispersion properties of modes 
localized at a single interface or in a slab of LH material 
[6-10]. 

In this paper, we present a comprehensive study of the 
properties of both linear and nonlinear surface waves at 
the interface between semi-infinite materials of two types, 
left- and right-handed ones, and demonstrate a number of 
unique features of surface waves in LH materials. In par- 
ticular, we show the existence of surface waves of both TE 
and TM polarizations, a specific feature of the RH/LH 
interfaces. We study in detail TE-polarized nonlinear 
surface waves and suggest an efficient way for engineering 
the group velocity of surface waves using the nonlinear- 
ity of the media. The dispersion broadening of the pulse 
can be compensated by the nonlinearity, thus leading to 
the formation of surface-polariton solitons at the RH/LH 
interfaces with a distinctive vortex-like structure of the 
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energy flow. We must note here, that the presented study 
is based on the effective medium approximation, which 
treats the LH materials as homogeneous and isotropic. It 
can be applied to the manufactured metamaterials, which 
possess negative dielectric permittivity and negative per- 
meability in the microwave frequency range, when the 
characteristic scale of the variation of the electromagnetic 
field (e.g. a field decay length and a wavelength of radia- 
tion) is much higher, than the period of the mctamaterial. 
The possibility of preparing isotropic LH materials was 
studied in [11], where the isotropy of the composite in 2D 
has been shown. To obtain a negative-refraction material 
in optics it is suggested that metallic nanowires are used 
Ref. [12]. Also, we note that losses is an intrinsic feature 
of LHM. However, the study of the effect of losses on the 
guided waves is not in the scope of the present paper. 

The paper is organized as follows. In Sec. II we study 
the properties of surface waves in the linear regime. We 
consider the most general case of an interface between lin- 
ear RH and LH media, and present a classification of TE- 
and TM-polarized surface waves localized at the inter- 
face. Section III is devoted to the study of the structure 
and general properties of nonlinear surface waves. We de- 
scribe the intensity-dependent properties of TE-polarized 
surface waves in three possible cases. In the first case, we 
assume that both the LH and RH media are nonlinear. 
In the second case, the LH medium is nonlinear, but the 
RH medium is assumed to be linear. In the third case, 
the RH medium is considered to be nonlinear while the 
LH medium remains linear. In all these cases, we take 
the nonlinear medium have an intensity-dependent Kerr- 
like dielectric permittivity. In Section IV, we study the 
frequency dispersion of nonlinear surface waves. In par- 
ticular, we demonstrate that the group velocity of surface 
waves can effectively be engineered by using the intensity- 
dependent dispersion. A detailed analysis is carried out 
for the example of nonlinear RH and linear LH media. 
Finally, in section V we describe the properties of non- 
linear localized modes propagating along the interface, 
and predict the existence of surface polariton solitons. 



II. LINEAR SURFACE WAVES 
A. Model 

Linear surface waves are known to exist, under certain 
special conditions, at an interface separating two differ- 
ent isotropic dielectric media. In particular, the existence 
of TM-polarized surface waves requires that the dielec- 
tric constants of two dielectric materials separated by 
an interface have different signs, whilst for TE-polarized 
waves the magnetic permeability of the materials should 
be of different signs (see, e.g. Ref. [13, 14] and references 
therein). Materials with negative e are readily available 
(e.g., metals excited below a critical frequency), whilst 
materials with negative /i were not known until recently. 
This explains why only TM-polarized surface waves have 



been of interest over the last few decades. 

In this paper, we consider an interface between the 
RH (medium 1) and LH (medium 2) semi-infinite me- 
dia, as shown in the inset of Fig. 1. The propagation of 
monochromatic waves with the frequency u is governed 
by the scalar wave equation, which for the case of the TE 
waves is written for the ^/-component of the electric field, 



dz 2 



+ 



dx 2 



-e(x)fj,{x) 



1 dfi(x) d 



fj,(x) dx dx 



Ey=0. 



(1) 

In the case of the TM waves, the scalar wave equations 
is written for the y-component of magnetic field, 



d 2 d 2 uo 2 . . . . 
^2 + ^2 + 



1 de{x) d 



e(x) dx dx 



H„ = 0. 



(2) 

In Eqs. (1) and (2), the functions e(x) and /j,(x) are di- 
electric permittivity and magnetic permeability in a bulk 
medium, respectively; ui is the angular wave frequency, 
and c is the speed of light in vacuum. The nonzero com- 
ponents of the magnetic field and of the electric field are 
found from the Maxwell's equations, i.e. for TE waves 



ic dE y 
u)fi dx 



and for TM waves 



E, 



ic 



dH„ 



dx ' 



ic dE y 
ojfi dz 

ic dH y 
ive dz 



(3) 



(4) 



respectively. Solutions of Eqs. (1) and (2) in each linear 
medium for localized modes, i.e. those propagating along 
the interface and decaying in transverse direction, have 
the form 



(E y ,H y )=A e ihz - K ^W, 



(5) 



where A is the wave amplitude at the interface, h is a 
propagation constant, 



«1,2 = 



h 2 - £l,2Ml,2 



1/2 



is the transverse wave number which characterizes the 
inverse decay length of the surface wave in the corre- 
sponding medium. 

It follows from Eqs. (1) and (2) that the tangen- 
tial components of the electric and magnetic fields 
change continuously at the interface between two media. 
These conditions give the dispersion relations for surface 
waves [6], 



^ + ^ = 0, 
Mi M2 



and 



Kl K2 

ei e 2 



0, 



(6) 



(7) 



for the cases of the TE- and TM-polarized surface waves, 
respectively. 
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FIG. 1: Existence regions of surface waves on the parameter 
plane (X, Y), where X — | £2 1 /ei and Y = | jl*2 | //xi - The inset 
shows the problem geometry. 




FIG. 2: Dispersion curves of the TE-polarized surface waves, 
for different values of ei , shown for the normalized values u> = 
cj/uip are h — hc/uj v . Dotted curves marks the dependence 
h = uj v /62H2- Dashed line is the critical frequency uii. 



B. Properties of surface waves 

For the analysis presented below, it is convenient to 
rewrite the dispersion relations (6) and (7) in the follow- 
ing form, 



and 



h 2 



h 2 



eiMi 



eiMi 



(?) 



Y(Y - X) 

(y 2 -i) : 



(8) 



2 X(X - Y) 
(X 2 -l) ' 



respectively, where we introduced the dimensionless nor- 
malized ratios X — |e 2 |/ei and Y — |/i2 which char- 
acterize the relative properties of the media creating the 
interface. The existence regions for surface waves can 
be determined from the condition of surface wave local- 
ization, i.e. when the transverse wave numbers are 
real, h > max {we\ji\/c, we2/i2/c}. Existence regions for 
both polarizations of surface wave are presented on the 
parameter plane (X, Y) in Fig. 1. Along with the polar- 
ization, we determine the type of the wave as forward or 
backward, as discussed below in Sec. II C. We note that 
there exist no regions where both TE- and TM-polarized 
waves co-exist simultaneously, but both types of surface 
wave can be supported by the same interface for different 
parameters, e.g. for different frequencies. 

One of the distinctive properties of the LH materials 
which has been demonstrated experimentally is their spe- 
cific frequency dispersion. To study the dispersion of the 
corresponding surface waves, it is necessary to select a 
particular form of the frequency dependence of the di- 
electric permittivity and magnetic permeability of the 
LH medium. A negative dielectric permittivity is se- 
lected in the form of the commonly used function for 
plasmon investigations [14] and a negative permeability 
is constructed in an analogous form (see, e.g. Ref. [6]), 
i.e. 



where losses are neglected, and the values of the parame- 
ters uip, uj r , and F are chosen to fit approximately to the 
experimental data [2]: lu p /2it = 10 GHz, w r /27r = 4 GHz, 
and F = 0.56. For this set of parameters, the region in 
which permittivity and permeability are simultaneously 
negative is from 4 GHz to 6 GHz. 

The dispersion curves of the TE-polarized surface wave 
(or surface polariton) calculated with the help of Eq. (10) 
are depicted in Fig. 2 on the plane of the normalized pa- 
(9) rameters ui = uj/lo p and h = hc/uj p . We note that the 
structure of the dispersion curves for surface waves de- 
pends on the relation between the values of the dielec- 
tric permittivities of the two media at the characteristic 
frequency, wi, at which the absolute values of magnetic 
permeabilities of two media coincide, fxi = \/j,2(lui)\. The 
corresponding curve in Fig. 2 is monotonically decreas- 
ing for ei > e2(wi)|, but it is monotonically increasing 
otherwise, i.e. for ei < |e 2 (o>i)|. Only the first case was 
identified in the previous analysis reported in Ref. [6]. 
The change of the slope of the curve (the slope of the 
dispersion curve represents the group velocity) with the 
variation of the dielectric permittivity of the RH medium 
can be used for group velocity engineering, which we dis- 
cuss in Sec. IV. 



e 2 (w) = 1 - 



IV. 



/j, 2 {oj) = 1 - 



Flu 2 



UJi. 



(10) 



The critical value of dielectric permittivity |e2(wi)| for 
the case of a nonmagnetic RH medium (fi\ — 1) is found 
from the dispersion relations (6) and (10), and it has the 
form 



e c = |e 2 (wi) = 1 



F 



- 1 



(11) 



For the parameters specified above, this critical value is 
e c = 3.5. 

The change of the dispersion curve from monotonically 
increasing to monotonically decreasing, shown in Fig. 2, 
is connected with a change in the direction of the total 
power flow in the wave, as discussed below. 
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C. Energy flow near the interface 

The energy flow is described by the Poynting vector, 
which defines the energy density flux averaged over the 
period T — 2it/uj, and can be written in the form 



S = —Re [E x H* 



(12) 



where E, H are the complex envelopes of the electric field 
and magnetic field of a surface wave, respectively, and the 
asterisk stands for the complex conjugation. 

A uniform surface wave propagating along the inter- 
face has only one non-zero component of the averaged 
Poynting vector, |S| = \S Z \. The energy flux in the RH 
and LH media is an integral of the Poynting vector over 
the corresponding semi-infinite spatial region, 



Pi 



f 



S z dz = 

Kl 



f 

Jo 



Bh 
dz = 

«2 



1/fjLi; for TE, 
1/ei; for TM, 



l//i 2 ; for TE, 
l/e 2 ; for TM, 



(13) 



(14) 



where the constant B = c 2 Aq / IQttui . We note that the 
electromagnetic energy flow is in opposite directions at 
either side of the interface, as was also predicted in [13] 
for TM polarized waves. The total energy flux in the 
forward z-direction is defined as the sum, P = Pi + P 2 , 
and it is found as 



P = (1 - XY) 



Bhu?e x[ i x J (1 + Y 2 ) j [y(/xiKi - // 2 k 2 )] , 

KlK2 c 2 { (l + X 2 ) /[X(e lKl -€2K 2 )], 

(15) 

for the TE and TM waves, respectively. The total energy 
flux is positive for XY < 1, and negative for XY > 1. 
The surface waves are forward or backward, respectively. 
The corresponding types of surface waves determined 
from this analysis are labelled in Fig. 1. 



III. NONLINEAR SURFACE WAVES 
A. Nonlinear LH/RH interface 

Nonlinear surface waves at an interface separating two 
conventional dielectric media have been analyzed exten- 
sively for several decades starting from the pioneering 
paper [15]. In brief, one of the major findings of those 
studies is that the TE-polarized surface waves can exist 
at the interface separating two RH media provided that 
at least one of these is nonlinear, but that no surface 
waves exist in the linear limit. 

In this section, we study TE-polarized nonlinear sur- 
face waves assuming that both media arc nonlinear, i.e. 
they display a Kerr-type nonlinearity in their dielectric 
properties, namely 



where the first term characterizes the linear properties, 
i.e. those in the limit of vanishing wave amplitude. 

First, we should mention that the recent systematic 
study of nonlinear properties of metallic composites [16] 
suggested the possibility of hysteresis-type nonlinear ef- 
fects in a structure consisting of arrays of split-ring res- 
onators (SRRs) and wires embedded in a nonlinear di- 
electric medium. Such effects can also be caused by a 
nonlinear dielectric material placed in the slits of the 
SRRs, which results in an intensity-dependent capaci- 
tance of the slit. These hysteresis effects can be avoided 
if the structure is filled by a nonlinear dielectric material 
except in the SRR slits. In what follows, we consider 
such composite structures for which the nonlinear prop- 
erties can be characterized by Eq. (16) valid far from the 
resonances. 

For a conventional (or right-handed) dielectric 
medium, positive an corresponds to a self-focusing non- 
linear material, whilst negative ot\ characterizes defocus- 
ing effects in the beam propagation. However, this clas- 
sification becomes reversed in the case of LH materials 
and, for example, a self-focusing LH medium corresponds 
to negative a 2 . Indeed, taking into account relation (16), 
we rewrite Eq. (1) for the case of the TE-polarized wave 
in nonlinear media as follows, 



d 2 E d 2 E 
dz 2 dx 2 



(en + ^a\E\ 2 ) E = 0. (17) 



According to Eq. (17), the sign of the product [ia de- 
termines the type of nonlinear self-action effects which 
occur. Therefore, in a LH medium with negative \i 2 all 
nonlinear effects are opposite to those in RH media with 
positive Hi, for the same a. Below, we assume for def- 
initcness that both LH and RH materials possess self- 
focusing properties, i.e. ct\ > and a 2 < 0. 

We look for the stationary solutions of Eq. (17) in the 
form Ei, 2 (x, z) = ^>i t2 (x) exp (ihz). Then, the profiles of 
the spatially localized wave envelopes ^i y2 (x) are found 
as [17] 



*i,2(z) = yii,2\] 2/011,21*1,2 J sechfai^ar-xi^)], (18) 

where 771.2 = Ki t2 c/uj arc normalized transverse wave 
numbers, xi, 2 are centers of the sech-functions which 
should be chosen to satisfy the continuity of the tan- 
gential components of the electric and magnetic fields at 
the interface. These conditions can be presented in the 
form of two transcendental equations, 



tanh 2 (7/1X1) 



1 



oi/zi?t|\ 
a 2 ii 2 ri 2 ) 



V OL2H1) 



and 



f]2 Vl 

— tanh (r] 2 x 2 ) = — tanh (771X1). 
M2 Mi 



(19) 



(20) 



: 1,2 



ei.2 + OiolEl 



(16) 



Note, that if the parameters {x\,x 2 ) correspond to one of 
the solutions of the equation, then (— Xi, —x 2 ) gives an- 
other solution. Two waves described by these solutions 
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have the same wave number, but they correspond to dif- 
ferent transverse structures of the surface wave, one of 
which has a maximum of the intensity at the interface, 
and the other one that has two humps shifted into the 
media. 

In a degenerate case, when (ai, /_ti, ei) = 
{\ct2 1 , M2 1 ? 1^2 1 ) : the system has an infinite number 
of solutions. Indeed, any (a;, — x) pair will describe a 
stationary solution for the surface wave. Such waves 
have a zero total energy flux because of the symmetry of 
the solution. 

The energy flow in this wave can be written in the form 



P = Po7 



aiMi 



M2 Ml 



"2^2 



tanh (771X1) 



«iMi, 

where Pq = c 2 /Airto 0*2^2, and 7 = hc/u is the normal- 
ized wave number. We now consider the surface waves 
in the non-degenerate case when only ct\ — |a 2 |. The 
dependence of the normalized energy flux P/Pq on the 
parameter 7 is shown in Fig. 3 for the cases when lin- 
ear waves are forward or backward, respectively. Cor- 
responding transverse wave structures are shown in the 
insets. 

The linear limit corresponds to the case P — ► when 
Xi — > +00 and X2 — > —00. Moving along the curves, the 
centers of the sech-functions move toward the interface 
and at the point with dP/drj = 00, x\ — X2 = 0, and 
from that point along the dashed line x\ — ► — exo, X2 —> 
+00, thus revealing the two-humped transverse structure 
of the surface wave. Note that the forward (backward) 
wave in the linear case remains forward (backward) in the 
nonlinear case, i.e. the type of the mode is determined 
by the linear parameters of the system and can be found 
using the diagram in Fig. 1. 



B. Nonlinear LH/linear RH interface 

Next, we consider surface waves propagating along an 
interface between linear RH and nonlinear LH media (see 
the inset in Fig. 4) having the negative nonlinear coeffi- 
cient «2 and, thus, displaying the self-focusing properties. 
The transverse structure of the stationary surface wave 
has the form: 

V(x) = { £oCxp ^^' x <0, 

\ {2/a 2 ^2) 1/2 m sech [r] 2 {x - x )} , x > 0, 

(22) 

where Eq and Xq are two parameters which should be de- 
termined from the continuity conditions at the interface 
for the tangential components of the electric and mag- 
netic fields, 



tanh {ri 2 x ) = M2»?i/(Mi»?2), 
Eq = (2/a 2 A*2) 1/2 ??2sech(?72Xo). 



(23) 



Analyzing the relations (23), we find that a surface wave 
always has the maximum of field intensity at the inter- 
face. This is in a sharp contrast to the nonlinear surface 




FIG. 3: Normalized energy flux vs. normalized wavenumbcr 
7 = hc/u for the nonlinear surface waves in two cases: (a) 
XY < 1, Y > 1, and (b) XY > 1, Y < 1. Solid curve corre- 
sponds to a one-humped structure, dashed - double-humped 
structure. The insets shows the structure of the surface waves 
at the points indicated by arrows. Dotted lines denote the lin- 
ear surface wave wavenumber. 



waves excited at the interface separating two RH media, 
when the electric field has maximum shifted into a self- 
focusing nonlinear medium [18]. 

The corresponding nonlinear dispersion relation of the 
surface waves is found in the form 



Mi M2 



1 - 



1/2 



0, 



(24) 



where A\ = i?o(ai/Ji/2) 1/ ' 2 is the normalized electric field 
amplitude at the interface. Equation (24) reduces to the 
linear dispersion relation (6) in the small-amplitude limit, 
i.e. when A\ — ► 0. 

The energy flux P associated with the nonlinear sur- 
face wave can be calculated in the form 



P = P0IV2 1 + 



Mi% 



where Pq is defined above. 



2 

M2 



1 



6 



Ph 



Y>1 



Y=l 



RH, 

lin&ar / 


l z LH 

. nonlinear 
\E(x) 


x o 





2.05 



Y<1 



2.1 2.15 
Y 




Linear limit 



2.2 2.25 



FIG. 4: Normalized energy flux vs normalized wavenumber 
7 = hc/co for the nonlinear surface waves at the nonlinear 
LH/Linear RH interface. Surface waves can be both forward 
(positive energy flux) and backward (negative energy flux). 
The inset shows the geometry of the problem. The solid line 
shows the transverse wave profile, dotted line shows the con- 
tinuation of the solution in nonlinear medium (22) to the lin- 
ear medium, dashed line indicates the position of the center 
xo of the sech - function. 




FIG. 5: Normalized energy flux vs normalized wavenum- 
ber 7 = hc/cu for the nonlinear surface waves at the lin- 
ear LH/nonlinear RH interface. Surface waves can be either 
forward (positive energy flux) or backward (negative energy 
flux). The inset shows the geometry of the problem. 



waves exist when Y > 1, reproducing the property of the 
corresponding linear waves. 



As an example, we consider the case XY < 1 for which, 
as we have shown above, only forward surface waves can 
exist at the interface between two linear media. However, 
nonlinear surface waves can be either forward or back- 
ward, as demonstrated in Fig. 4. For Y < 1, there exists 
no linear limit for the existence of the surface waves, 
while in other two regions the results for linear surface 
waves are recovered in the limit P — > 0. The point on 
the curve corresponding to P = in the case Y < 1 de- 
scribes the wave of finite amplitude in which the energy 
flows on either side of the interface are balanced. Such a 
wave does not exist in the linear limit. 



C. Linear LH/nonlinear RH interface 

Finally, we consider the case when the LH material 
is linear, while the RH medium is nonlinear. In such 
a geometry, the dispersion relation for the TE-polarized 
waves has the form 

=0, (25) 

where A2 = Eq yJcti^l'Z is the normalized amplitude of 
the electric field at the interface. 

The dependence of the normalized energy flux on the 
wave number of the surface wave is shown in Fig. 5 for 
XY > 1. In contrast to the linear waves, the nonlin- 
ear surface waves can be either forward or backward (see 
Fig. 1). In analogy with the case Y < 1 for the nonlin- 
ear LH/linear RH interface, it can be shown that there 
exists no small-amplitude limit for the nonlinear surface 
waves for Y > 1 . For XY < 1 only forward travelling 



IV. FREQUENCY DISPERSION OF 
NONLINEAR SURFACE WAVES 

We have demonstrated in Sec. II B that the frequency 
dispersion of surface waves depends on the dielectric per- 
mittivity of the RH medium (see Fig. 2). These results 
suggest that the dispersion type can be switched between 
normal and anomalous if the RH medium is nonlinear. 

To demonstrate this property, we study the properties 
of nonlinear surface waves near the critical point and se- 
lect ei = 3.4, in order to stay just below the critical value 
e c corresponding to the linear case. It should be men- 
tioned here that although the negative permeability of 
the LH composite material is necessary for the existence 
of TE surface waves in the present model, it has been 
shown [19] that TE surface waves do exist at the inter- 
face between a RH nonlinear plasma and a RH nonlinear 
dielectric medium with a constant zero-field permittivity, 
provided that the nonlinear parameter in the plasma ex- 
ceeds that in the dielectric medium. In that case, since 
(i = 1 in both media, it has also been shown [19] that 
the relation for wave intensity at the boundary and the 
frequency are independent of the wave number. In the 
present case, where the LH composite material has a per- 
meability not equal to unity, this result is no longer valid. 
In our problem, Eq. (25) provides a dependence between 
the three variables h, A 2 and u, so that the wave intensity 
at the boundary depends on both h and lo. 

As was shown for the case of linear surface waves, a 
change of the slope of the dispersion curve takes place 
at the critical value of dielectric permittivity of the RH 
medium (11). In the nonlinear case, the dielectric permit- 
tivity of the RH medium depends on the field intensity 
and, in particular, it exceeds the critical value when the 
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FIG. 6: Normalized frequency vs. normalized wave number 
for different values of the amplitude A 2 . Dashed line corre- 
sponds to the line uj\ in Fig. 2; dotted line is 7/2 = 0. 



wave amplitude A 2 becomes larger than the threshold 
value A 2c given by the equation 



A 2 



(26) 



For the media parameters corresponding to Fig. 2, equa- 
tion (26) gives: A 2c = 0.3162. Figure 6 shows the disper- 
sion curves for three different wave amplitudes. As was 
shown in Sec. II B, the frequency dispersion is normal for 
the wave amplitudes below the critical value (26), and it 
is negative, otherwise. One can also notice from Fig. 6 
that the existence region for surface waves depends on 
the wave amplitude A 2 . In Fig. 7 we present this region 
of wave existence on the plane of the wave amplitude and 
normalized frequency. The existence region of the back- 
ward surface waves below the critical value collapses at 
A 2c , and it expands in the region of the forward surface 
waves above the threshold. Note, that the forward and 
backward waves exist at different frequencies only. 

The existence regions of the surface wave can be ex- 
plained from the viewpoint of the physics of wave local- 
ization. The wave localization is determined by the nor- 
malized transverse wavenumbers 771,2, which define the 
inverse decay length of the surface wave in the corre- 
sponding medium. The conditions rji — and rj 2 = 
correspond to the delocalized waves, and determine the 
boundary of the existence region. Figure 8 shows the 
dependence of the normalized wave number on the wave 
amplitude for different frequencies. These curves rep- 
resent the horizontal cross-sections of the region of the 
wave existence shown in Fig. 7. The dashed line in Fig. 8 
shows the boundary of localization of the wave in the LH 
material. Comparing Fig. 8 and Fig. 7, we come to the 
conclusion that in Fig. 7 the upper boundary for the ex- 
istence of the forward waves (below A 2c ) and the lower 
boundary for the backward wave existence (above A 2c ) 
are determined by the wave localization in the LH mate- 
rial. 



FIG. 7: Existence region of nonlinear surface waves (shaded). 
The curves show the amplitude A2 vs. normalized frequency, 
for different values of the normalized wave number (marked 
at the curves) . All curves intersect at the critical point A2 = 
A 2c - 
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FIG. 8: Normalized wave number vs. normalized field ampli- 
tude A2, for different frequencies. Dashed line is 7/2 = 0. The 
curves on the left of A 2c meet the wave number axis at the 
values corresponding to the linear case. There exist no linear 
solutions to the right of the critical vertical line A2 = A 2c 



given by the result 



c 2 7 F 2 

--Hin , 



(27) 



and, using the boundary value technique [18], the power 
flow in the nonlinear half-space (RH medium) can be 
obtained in the following form 



Pi 



2 

c 7 



2 *UH 

2 



(28) 



We note here that, although the nonlinear coefficient evi 
appears as a factor in the denominator of Eq. (28), the 
reduction to the linear case (when a\ — ► 0) can be per- 
formed in a straightforward way by expanding Eq. (28) 
as for small ot\ as follows, 



Pi 



2 

c z 7 



16TTu>fj,ir]i 



El + 0{a\), 



(29) 



The power flow in the linear LH composite medium is 



which has the same form as Eq. (27). 
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FIG. 9: The absolute value of the ratio of the power flow in 
the nonlinear dielectric to that in the LH material. Above the 
critical value A c , the power flow in the nonlinear half-space 
dominates giving a forward travelling wave. As K2 approaches 
zero towards the left-hand side of the curves, the power in the 
left-handed medium dominates. 



The absolute value of the ratio of the power flow in 
the RH nonlinear medium to the power flow in the LH 
composite medium is depicted in Fig. 9. Above the crit- 
ical value A c , there exists a value of the wave number at 
which the power flow is positive. In this region, there ex- 
ist a forward travelling surface wave. Note that no matter 
what the value of the field intensity at the boundary is, 
there always exists some value of the wave number where 
the flow is negative and there exists a backward travelling 
surface wave. This can be seen by reference to Eq. (27). 
The presence of 772 in the denominator of Eq. (27) means 
that as h approaches a value that makes i] 2 = (i.e. 
a dispersion curve in Fig. 6 approaches the dashed line) 
the negative power flow in the LH material dominates the 
total power flow. Conversely, as 772 increases from zero, 
the negative power flow in the LH medium decreases so 
that, provided the intensity of the electric field at the 
boundary is high enough, the positive power flow in the 
nonlinear dielectric dominates. 



V. NONLINEAR PULSE PROPAGATION AND 
SURFACE- WAVE SOLITONS 

A. Envelope equation 

Propagation of pulses along the interface between RH 
and LH media is of a particular interest, since it was 
shown before [13] for TM modes that the energy fluxes 
are directed oppositely at cither side of the interface. 
Therefore, we can expect that the energy flow in a pulse 
with finite temporal and spatial dimension has a nontriv- 
ial form [10] and, in particular, it can be associated with 
a vortex-like structure of the energy flow. 

We analyze the structure of surface waves of both tem- 
poral and spatial finite extent that can exist in such a 
geometry. To obtain the equation describing the pulse 
propagation along the interface, we look for the structure 



of a broad electromagnetic pulse with carrier frequency 
loq described by an asymptotic multi-scale expansion with 
the main terms of the general form 



ifiQZ — iojQt 



[*o{x)A{Z,t)- 



i^i(x) 



*2(*,£,t) 



(30) 



where the field 'J? stands for the components (E y , H X ,H Z ) 
of a TE-polarized wave, the first term *o = 
(E y o, H x o, H z q) describes the structure of the mode at 
the carrier frequency ivq, 5*1 is the first-order term of the 
asymptotic series which can be found as ^1 = d^o/dh, 
and ^2 is the second-order term. Here A is the pulse 
envelope, h is the wave number corresponding to the 
carrier frequency u>o, £ = z — v g t is the pulse coordinate 
in the reference frame moving with the group velocity 
v g = duo/dh. Substituting Eq. (30) into Eq. (17) and us- 
ing the Frcdholm alternative theorem [20] , one can obtain 
the equation for the evolution of the field envelope 



8A 8d 2 A n , lA2lA n 
+ v 2 (h)\A 2 \A = 0, 



dt 2 de 



(31) 



where the coefficient 6 = d 2 u/dh 2 stands for the group- 
velocity dispersion (GVD) which determines the pulse 
broadening and can be calculated from the dispersion 
relations, u) 2 (h) = (dojNL/dA 2 ) \a=o is the effective non- 
linear coefficient calculated with the help of the nonlinear 
dispersion relation. The NLS equation has a solution in 
the form of a bright soliton localized at the interface, 
provided the GVD (6) has the opposite sign to the sign 
of the nonlinear coefficient (ui 2 ) (see, e.g., [21] and ref- 
erences therein). The existence of the surface polariton 
solitons has been predicted in a number of structures sup- 
porting nonlinear guided waves (see, e.g., Ref. [22] and 
references therein). 

The effective nonlinear coefficient for the case of an in- 
terface between the nonlinear LH medium and the linear 
RH medium, can be found using Eq. (23), 



u) 2 (h) 



duj 



4/ic 2 (e 2 /J2 — d/Ji) dh 



(32) 



The signs of the group velocity dw /dh and of the param- 
eter S can be determined from the Fig. 2. As a result, 
for any reasonable values of dielectric permittivity and 
magnetic permeability of the RH medium, there exists a 
range of frequencies for which lo 2 • 5 < 0, indicating the 
possibility of exciting surface polariton solitons. 

To study the energy flow in such a surface-polariton 
soliton, we use the asymptotic expansions (30) for the 
field components, and from Eq. (12) we obtain the energy 
flow structure described by their components 



S z = 



C h E 2 |^|2 



(33) 
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FIG. 10: A vortex-like distribution of the Poynting vector in a 
surface wave soliton propagating along the LH/RH interface. 



Re A 



,dA* 



(34) 



The structure of the Poynting vector in the surface-wave 
pulse is shown in Fig. (10), where it is clearly seen 
that the energy rotates in the localized region creating a 
vortex-type energy distribution in the wave. The differ- 
ence between the Poynting vector \S Z | integrated over the 
RH medium and that calculated for the LH medium de- 
termines the resulting group velocity of the surface wave 
packet. 

We note the distinctive vortex-like structure of the sur- 
face waves at the interface separating RH and LH media 
follows as a result of the opposite signs of the dielectric 
permittivity for the TM-polarized waves and the mag- 



netic permittivity for the TE-polarized waves. These 
conditions coincide with the conditions for the existence 
of the corresponding surface waves and, therefore, surface 
polaritons should always have such a distinctive vortex- 
like structure. 



VI. CONCLUSIONS 

We have presented a systematic study of both linear 
and nonlinear surface waves supported by an interface 
between a left-handed metamaterial and a conventional 
dielectric medium. For the linear regime, we have ex- 
tended some earlier theoretical results and analyzed dif- 
ferent types of surface waves and their existence regions. 
In particular, we have demonstrated that the structure 
of the energy flow in a spatially localized wave-packet of 
the surface waves propagating along the interface has a 
vortex-like structure. For the case of nonlinear surface 
waves, we have demonstrated that when only one of the 
media is nonlinear the maximum amplitude of the spa- 
tially localized wave does not shift out the interface. This 
result is in a sharp contrast to the case of an interface be- 
tween linear and nonlinear right-handed materials where 
the maximum of a surface wave is always shifted into a 
self-focusing nonlinear medium. We have demonstrated 
that the group velocity of nonlinear surface waves can be 
controlled by changing the intensity of the electromag- 
netic field and, in particular, the surface wave can be 
switched from the forward propagating one to the back- 
ward propagating one by varying the field intensity only. 
In addition, we have obtained the conditions for the ex- 
istence of the surface-polariton solitons at the metama- 
terial interface. 
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